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^vq , Abstract. The spectral functions for calcium and argon are constructed. It is verified that their 

predictions for the quasielastic electron-nucleus cross sections in the energy range ~1 GeV agree 
with the data. The argon spectral function is then used to obtain the quasielastic neutrino-nucleus 
. cross section. 
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The description of lepton-nucleus interactions in the 1 GeV energy region is based on 
the impulse approximation (IA). The most reliable implementation of the IA is that of 
the spectral function (SF). Exact calculations of the SF exist only for the lightest nuclei, 
whereas well founded approximations for carbon and oxygen [Ql 0]. In this note, we 
report results on modeling of the SFs for heavier nuclei like 2(jCa and jgAr. 

Our primary motivation was to provide a good description of medium-sized nuclei by 
their SFs and, as a result, to improve accuracy of the predicted neutrino-nucleus cross 
sections with respect to the Fermi gas (FG) model for neutrino energy around 1 GeV. 
We focus on quasielastic (QE) scattering and therefore verify obtained SFs in the region 
of energy transfer up to the QE peak. 

G\ ' 

Model description. Every SF is divided into two parts: the correlated and the mean 
field (MF) one. Analytical formula for the correlated part of the SF is given in Ref. dH]. 
The MF part is modeled on the simplifying assumption that the level widths does not 
depend on momentum. The obtained expression contains a contribution from each shell- 
r> \ model state a of energy E a occupied by N" nucleons: 
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/> MF (p,£) = i^fe)^ JXF„(£« -E) . (1) 

In the above, the MF part of the momentum distribution is denoted by nf^ and the 
number of nucleons of isospin t by N t . As the function F a , describing width of level a, 
we use the Gaussian distribution. For further details, see Q4£| . 

In the numerical computations, we use the momentum distributions from Ref. 
and recent parametrization of the electromagnetic form factors ML To evaluate the off- 
shell hadronic current matrix elements, de Forest's prescription [7] is adopted. In the 
electromagnetic case, a procedure to restore the current conservation is applied. 

Two separate final state interactions (FSI) effects are included: Pauli blocking and 
reinteractions of the struck nucleon described by means of the optical potential Jil]. 
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FIGURE 1. Differential cross section of the process 16 0(e, e') for fixed scattering angle. Results for our 
model (solid line) are compared to the Benhar SF [2] with the same FSI (dashed line) and to the Fermi 
gas model (dotted line). The experimental data are from The values of momentum transfer at the 

peaks are 638 MeV (for beam energy 1200 MeV), 441 MeV (for 730 MeV), and 325 MeV (for 537 MeV). 
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FIGURE 2. Differential cross section of 40 Ca(<?,<?') scattering at angle 45.5°. Calculations within our 
model (solid line) are compared to the results of Butkevich and Mikheyev UlTl (dashed line), the Fermi 
gas model (dotted line), and the experimental data [ 12]. The corresponding values of |q| at the peaks are 
602 MeV (for beam energy 841 MeV), 490 MeV (for 681 MeV), and 342 MeV (for 471 MeV). 



Logic of our work. The basic steps of our research were: 

• Identification of the region in the energy transfer and momentum transfer (co, |q|) 
plane, from which most of the neutrino QE cross section comes from. 

• Selection of the electron scattering data with similar 00 and |q|. 

• Comparison of the results for oxygen obtained within our model to the available 
electron scattering data [EIQli and the Benhar SF OL0]. 

• Verification of accuracy for ^Cafe, e') and comparison with the results of Ref. ifllll . 

These two tests showed satisfactory agreement of our approach with the data, see Figs.Q] 
and El Therefore the next step was taken: Application of the model to argon. The 
left panel of Fig. [3] presents comparison with the few known experimental points for 
jgAr(e,e') scattering fl3ll . Since the accuracy is comparable to the case of oxygen [0], 
we proceed to our final goal and apply the new SF to neutrino scattering. The typical 
plot is presented in the right panel of Fig. [3l 
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FIGURE 3. Differential cross sections of argon. Comparison of our approach (solid line) to the FG 
model (dotted line). Left panel: The process jgAr(e,e'). Experimental points are from j 1311 . The value of 
momentum transfer at the peak is 371 MeV. Right panel: The process igAr(Vjj, 



Discussion. Only the QE scattering was considered. Thus for higher values of ft), 
some of the cross section is lacking. The oxygen nucleus provides the opportunity 
to compare our approximation with the more systematic Benhar SF, see Fig. CD The 
simplifying assumptions are responsible for slightly different shape of the QE peak 
obtained within our model, but the discrepancy is of the size of the error bars, and the 
achieved accuracy is good. Figure [2] shows that for heavier calcium target, the agreement 
with data is better than for oxygen. 

We conclude that the way we model SFs works fine in the region of ft) and |q| relevant 
to ~ 1 GeV neutrino QE scattering. By adopting the standard extrapolation, we expect 
that the model tested for electrons describes neutrino scattering equally well. It allows us 
to presume that the accuracy of the neutrino cross section plot shown in the right panel 
of Fig. [3]is similar to those presented for electron scattering. 

In the case of all three targets, larger deviations from the electron scattering data occur 
for lower beam energies. This can be understood as an indication that when momentum 
transfer at the QE peak is lower than ~400 MeV, the IA is no longer reliable. 
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